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This talk presents results from the QCDSF-UKQCD collaboration for moments of leading twist generalized 
parton distributions in two-flavor lattice QCD based on 0(a) improved Wilson Fermions. We study helicity 
independent and helicity flip GPDs with a focus on densities of quarks in the transverse plane. 


1. INTRODUCTION 


and the tensor GPDs are defined by [5,6] 


Generalized parton distributions (GPDs) [1,2] 
enable us to study many fundamental aspects of 
the intrinsic hadron structure in the framework of 
QCD. They seamlessly link the at first sight com¬ 
pletely different concepts of parton distributions 
and form factors. Furthermore, they allow us to 
define and calculate orbital angular momentum 
contributions of quarks and gluons to the nucleon 
spin [3]. Particularly important for this study is 
the interpretation of GPDs as probability densi¬ 
ties in impact parameter space [4]. The GPDs H 
and E of quarks in the nucleon are defined by the 
following off-forward nucleon matrix elements 

(P',A'| |P,A) = 

+ ( 1 ) 
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{P\A'\J^E^^q(-^r^n^a>^''j5l^qyyP,A) = 


u{P',A')nfj,\a>^‘'j5{HT{x,C,t) - -^HTix,^,t) 


2m 


+ 


2m? 


ETix,tt) 


ETix,^,t)\u{P,A) , 


( 2 ) 


where A = P is the mo¬ 

mentum transfer with t = A^, P = (P' -|- P)/2, 
and ^ = — n • A/2 defines the longitudinal mo¬ 
mentum transfer with the light-like vector n. The 
Wilson line ensuring gauge invariance of the bilo¬ 
cal operator is denoted by U. Here and in the fol¬ 
lowing we do not show explicitly the dependence 
of the GPDs on the resolution scale Q^. 

A prominent feature of GPDs is that they re¬ 
produce the well-known parton distributions in 
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the forward limit, A = 0. The relations to the un¬ 
polarized q{x) and transversity 6q(x) parton dis¬ 
tributions are given by H{x,0,0) = q{x) = fi{x) 
and Ht{x, 0,0) = Sq{x) = hi{x), respectively. 
On the other hand, integrating the GPDs over x 
gives the form factors: 


J dxH{x,i,t) = Fi(t), 

J dxHT{x,^,t) = grit). 


( 3 ) 


Another property of GPDs important for our in¬ 
vestigations below is their interpretation as den¬ 
sities in the transverse plane for ^ = 0 [4]. To give 
an example, it has been shown that the impact 
parameter dependent quark distribution for the 
quark GPD Hg, 

qix,b±) = J =-Ai), 

( 4 ) 


has the interpretation of a probability density 
for unpolarized quarks of flavor q with longitudi¬ 
nal momentum fraction x and transverse position 
b± = {bx, by) relative to the center of momentum 
in a nucleon. 


2. LATTICE CALCULATION OF 
GENERALIZED FORM FACTORS 

In order to facilitate the computation of the 
GPDs in lattice QGD, we first transform the LHS 
of Eq. H1I2II to Mellin space by forming the in¬ 
tegral dxx"~^ ■ ■ ■. This gives off-forward nu¬ 
cleon matrix elements of towers of local operators 

= q(0)a>^^''-f5iD^^ ... , 

( 5 ) 

where D = ^{D — D) and {• • •} indicates sym- 
metrization of indices and subtraction of traces. 
Matrix elements of the operators in Eq. are 
parametrized in terms of generalized form fac¬ 
tors (GFFs) Ani, Bni and Cno for the vector, and 


Ar-ni, BTni, Axm and Bxni for the tensor case. 
For the vector operators, the parametrizations for 
the lowest three moments read 

(P'A'|9(0)7'^g(0)|PA) = 

r I 

u{P',A')h>^AMt) + j^r(P, A) , 

( 6 ) 

(P'A'I 9(0)7foiB'"i^g(0) |PA) = 

uiP', A') I> A 20 it) + B 20 it) 

+lAfoA'^i>C2o(t)U(P,A) , (7) 

m j 


and 


(P'A'I 9(0)7fozB'"GB'"">g(0) |PA) = 

Il(P',A')|p^''p''^7'^=>A3o(t) 

^^ AfoA^i7'^=>A32(t) 
2m 

-bAfoA^^ •' Bs2it)'>uiP,A) . (8) 


For the n = 1 helicity flip operator, we have [5,7] 
(P'A'I 9(0)a'^'^759(0)|PA) = 
m(P', A')|cr'^'^75 




2m 


A['^Cr^l“ 75 Ao 


2m^ 


ATioit)>uiP,A) . 


( 9 ) 


Parametrizations for higher moments n > 1 can 
be found in [7,8]. 

Lattice counterparts of the continuum 
Minkowski space-time matrix elements in Eqs. 

0 are nucleon two- and three-point correlation 
functions within a discretized Euclidean space- 
time framework. They are defined by 


C2P*(r,P) = ^ (^) (A^fc(T,P)A,(r.rc,P)) , 

jk ^ ''jk 


( 10 ) 
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and 


J2^jk{NkiTsnl.,P')0^^-^-^{T)N,iT,,c,P)}, 

j,k 

( 11 ) 

where F is a (spin-)projection matrix and the 
operators N and N create and destroy states 
with the quantum numbers of the nucleon, re¬ 
spectively. Inserting complete sets of states and 
using the time evolution operator, we can rewrite 
Eq. itTTl) . 

C'oPt^'-''"-i(r,P',P) = 

{Z{P)Z{P')) E(P)T-E(P')(T.„^-r) 

AE{P')E{P) 

X |P,A) 

A,A' 

x?I(P,A)frt(P',A') + ... , (12) 

where we set Tgrc = 0. The two-point function 
can be written as 


(^2pt(T,p) ^ 

+ (13) 

Here and below, m denotes the mass of the nu¬ 
cleon ground state. Excited states with energies 
E' > E{P),E{P') are represented by the ellip¬ 
sis in Eq. m and m- They are exponentially 
suppressed as long as t l/A’^ Tgnk — 1/E' . 

Using the explicit parametrizations from Eqs. 

EJ transformed to Euclidean space together with 
Eq. flUl . we can sum over polarizations to obtain 


c: 


P',P) = 


AE{P')E{P) 


f{ip' - m){ a^^ -'"”-M„o(f) 


xTr 


+&''i ''^"”iPno(i) H- (*/’-w) 


(14) 


where a,b... are the Euclidean versions of the 
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Figure 1. Effective dipole masses for the iso¬ 
vector form factor Aio = Pi for all datasets. 


corresponding prefactors in Eqs. (EEJ. The 
Dirac-trace in Eq. 03 is evaluated explicitly, 
while the normalization factor and the exponen¬ 
tials are cancelled by constructing a ratio R of 
two- and three-point functions, 

Ro{t,P\P) = 

■C2pt(^^p/)C2pt(r,„k,P')(^"P*(7-snk-r,P)1 ^ 

"" C2pt(r,P)C2pt(^^„^^P)C'2pt(r,„k-r,P') ' 

(15) 

The ratio R is averaged over the plateau region 
and then equated with the corresponding sum 
of GFFs times P- and P'-dependent calculable 
pre-factors. For a given moment n, this is done 
simultaneously for all contributing index combi¬ 
nations (/ii.. .p,„_i) and all discrete lattice mo¬ 
menta P, P' corresponding to the same value of 
t = (P' — P)^. Following this procedure, we 
find in general an overdetermined set of equa- 
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Figure 2. Slopes of the generalized form factors 
Aio{t) (lowest curve), A 2 o{t) and A 3 o{t) (highest 
curve). 


Figure 3. Decrease of the charge radius with in¬ 
creasing n. 

disconnected contributions. We expect, however, 
that they are small for the tensor GFFs [13]. 


tions from which we extract the GFFs [9,10]. On 
the lattice the space-time symmetry is reduced 
to the hypercubic group H{4), and the lattice 
operators have to be chosen such that they be¬ 
long to irreducible multiplets under H{A). For 
our analysis of the n = 1,2,3-vector operators, 
we employ all operators belonging to the corre¬ 
sponding 4- , 9- and 12-dimensional spaces de¬ 
scribed in [11], and for the n = 1 tensor opera¬ 
tors we utilize the 6-dimensional multiplet as dis¬ 
cussed in [6]. The statistical error on the GFFs 
is obtained from a jackknife analysis. Our results 
have been non-perturbatively renormalized [12] 
and transformed to the MS scheme at a scale of 4 
GeV^. The lattice results to be discussed below 
have been obtained from simulations with n/ = 2 
flavors of dynamical non-perturbatively 0(a) im¬ 
proved Wilson fermions and Wilson glue. There 
are 12 datasets available consisting of four dif¬ 
ferent couplings (3 = 5.20, 5.25, 5.29, 5.40 with 
three different n = Ksea values per /3. The pion 
masses of our calculation vary from 550 to 1000 
MeV, and the lattice spacings and spatial volumes 
vary between 0.07-0.11 fm and (1.4-2.0 fm)^ re¬ 
spectively. Our simulations are based on three 
sink momenta, po = (0,0,0), pi = (p, 0,0), 
P 2 = (0,p, 0) {p = 27r/Ls) and three projection 
matrices, funpoi = ^( 1 -^ 74 ), Fi = ^( 1 + 74 ) *757i: 
F 2 = ^(1 + 74 )^ 7572 - Our present calculation 
does not include the computationally demanding 


3. LATTICE RESULTS FOR THE 
LOWEST MOMENTS OF QUARK 
DENSITIES IN THE NUCLEON 


In a finite volume, momenta are discretized 
and we obtain the GFFs only for a limited num¬ 
ber of different values of the momentum transfer 
squared t. In our simulation, there are in gen¬ 
eral 16 t-values available per dataset in a range 
of 0 < t < 4 GeV^. In order to get a represen¬ 
tation of the GFFs for continuous values of t, we 
parametrize them using a p-pole ansatz 


m = 


m 


(16) 


where the values of E( 0 ), irip and p for the in¬ 
dividual GFFs are fixed by a fit to the lattice 
results. The ansatz in Eq. ca is then Fourier 
transformed in order to get the GFFs in coordi¬ 
nate space. For the vector GFFs Aio, A 20 and 
A 30 , we use a dipole ansatz, p = 2. More details 
of the p-pole parametrization and some numerical 
results for the parameters can be found in [16-18]. 
Let us note that the values for the power p we are 
using for the different GFFs lead to a regular be¬ 
havior of the quark densities in the limit b± 0 , 
as discussed in [19]. In order to check our dipole 
ansatz, we invert Ea. (ll6ll and show in Fig.(Q the 
effective dipole-mass squared = m.p= 2(0 
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Figure 4. Spin independent densities of up-quarks in the nucleon for n = 1, 2,3 (from left to right). 


as a function of the momentum transfer squared 
t for the x4io(t) = Fi{t) form factor. The hori¬ 
zontal lines indicate the results for from fits 
to the individual datasets using Eq. The 

plot shows that the dipole represents an accurate 
description of the form factor for a large range of 
the momentum transfer. 

A momentum transfer squared of the order of 
a few GeV^ or more has in general to be redis¬ 
tributed between the partons in order to prevent 
a break-up of the nucleon. According to the Feyn¬ 
man mechanism, parton configurations for which 
the active quark has a large longitudinal momen¬ 
tum fraction may still form a bound state after 
the scattering without a redistribution of the mo¬ 
mentum transfer. This simple picture of the scat¬ 
tering mechanism is confirmed by our results in 
Fig.@, where we compare the slope in t of the 
generalized form factors A„o for the lowest three 
moments, which have been normalized to unity at 
t = 0 GeV^. To see this most clearly, let us first 
note that higher moments n correspond to larger 
average momentum fractions x" « (a:")/(a;”“^) 
[20]. Fig.@ shows that already at a momentum 
transfer squared of t ~ — 2 GeV^ the (normalized) 
form factor A 30 is a factor of two larger than the 
form factor Aio. This indicates that a momentum 
transfer of this magnitude to the nucleon can be 
better absorbed by partons with large momentum 
fraction x, thereby lowering the probability of a 
break-up of the bound state. In the extreme case 
n —*■ 00 , we expect the generalized form factors to 


be independent of t, A„^oo.o(f)Mn^cx3.o(0) = 1- 
The flattening of the slope for the A„o GFFs in 
lattice QCD has for the hrst time been observed 
in [14]. 

In impact parameter space, the increase of the 
dipole masses translates into a decrease of the 
charge radius for increasing average momentum 
fraction. Fig. m shows that the charge radius of 
the nucleon in the transverse plane drops by a 
factor of about two, going from n = 1, which 
corresponds to an average momentum fraction of 
« 0.25 to n = 3, where « 0.4. Since 
the impact parameter measures the distance of 
the active quark to the center of momentum R ±, 
we expect the distribution of the quarks in the 
transverse plane to peak around i?j_ = 0j_ when 

X ^ 1 [21]. 

The x-moments of the impact parameter de¬ 
pendent quark distribution in Eq. 0 are given by 
the GFFs in impact parameter space, 

J x'^~^q{x,b±) = Ano{b±). (17) 

Fourier-transforming the p-pole ansatz for the 
GFFs in Ea. dltill to impact parameter space, we 
can calculate the distribution of quarks in the 
transverse plane for the lowest three moments 
we consider in this work. The result is shown 
in Fig. 0 , where we plot densities of up-quarks 
in Eq. 03 for n = 1 , 2 , 3 . To be precise, let us 
note that the plots do not exactly show proba¬ 
bility densities because the moments x"“^ al- 
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Figure 5. Densities of up-quarks in the nucleon. The nucleon and quark spins are oriented in the 
transverse plane as indicated, where the inner arrow represent the quark and the outer arrow the nucleon 
spin. A missing arrow represents the unpolarized case. 


ways correspond to differences/sums of quark and 
anti-quark densities. The quark distributions in 
Fig.® clearly get narrower from the left (n = 1) 
to the right {n = 3) plot. 

In the following, we will discuss the recently 
observed correlation of transverse quark spin and 
impact parameter which shows up in transverse 
spin densities of quarks in the nucleon [19]. It 
turns out that these correlations in the trans¬ 
verse plane are governed by quark helicity flip (or 
tensor) GPDs. The lowest moment of the quark 
transverse spin density reads [19] 

(p+,i?j_ = 0,5j_|/5r'(^±)|^"+,i?± = 0,5^) = 

+ sl(2516i - , (18) 

where the lowest moment of the transverse spin 
density operator is given by 

= ^q{b±)['y~^ - si ter 75 ] g(&j.) (19) 

which projects out quarks with transverse po¬ 
larization s_L = (cosx', siny'). The transversity 


states 


P+,Pj. = 0,^i) = ^(|P+,Pi 
+e*^|P+,Pi = 0,A=-)) 


0, A — -|-) 


( 20 ) 


describe a nucleon with longitudinal momentum 
P+ = (P° -I- P^)/\/2 which is localized in the 
transverse plane at R± = 0 and has trans¬ 
verse spin S± = (cosy,sinx). The impact pa¬ 
rameter dependent GFFs in Eq. m are the 
Fourier-transforms of the momentum space GFFs 
at ^ = 0, as in Eq. 0 - The derivatives in 
Eq. llTKIl are defined by f'{b±) = 81 , 2 ^f{b±) and 
Ab^f{b±) = 4 db 2 ^{b'^^db 2 jf{b±). In Fig.® we 
show preliminary results for the lowest moment of 
transverse spin densities of quarks in the nucleon 
for up quarks. The densities are strictly positive 
for all b±, indicating that the contributions from 
anti-quarks are small. On the LHS of Fig. ®, 
we show the density of unpolarized quarks in a 
transversely polarized nucleon. The distortion of 
the density in -t-y-direction coming from the term 
cx P[o(^J-) already been discussed in [21]. 
In [22], a formalism has been developed which 
predicts that such a distortion leads to single spin 
asymmetries and in particular to a non-zero neg¬ 
ative Sivers function [23] /j^ < 0 for up-quarks. 
First results from experiment seem to confirm 
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this [24], The density plot in the center of Fig.0 
shows that the term oc in Ea. ltT^ also 

leads to a strong transverse deformation orthog¬ 
onal to the transverse quark spin for an unpolar¬ 
ized nucleon. According to [15], this distortion in 
-|-j/-direction^ due to the transversely polarized 
quarks may correspond to a non-zero, negative 
Boer-Mulders function [25] < 0 for up-quarks. 

For the plot on the RHS in Fig. ©, the distortions 
due to transverse quark and nucleon spin add up, 
while there is practically no influence visible from 
the quadrupole-term oc s[^(26[^5^ — b']_S’’^)Si_ in 
Eq. itT^ . 

4. CONCLUSIONS AND OUTLOOK 

We presented results for the lowest moments of 
helicity independent and helicity flip generalized 
parton distributions in ny = 2 lattice QCD. Al¬ 
ready at a momentum transfer of a few GeV^, we 
find that the vector coupling to the nucleon is no¬ 
ticeably stronger for quarks with larger momen¬ 
tum fraction. Equivalently, the charge radius of 
the nucleon decreases visibly with increasing av¬ 
erage momentum fraction. For transverse polar¬ 
izations, we observe strong distortions of the spin 
densities orthogonal to the direction of the quark 
and nucleon spin. These distortions could give 
rise to non-vanishing Sivers and Boer-Mulders 
functions through final state interactions as ar¬ 
gued by Burkardt [15,22]. 

We plan to extend the analysis of moments 
of generalized parton distributions, including a 
study of the lowest two moments of transverse 
spin densities for up- and down-quarks, and to 
investigate improved positivity bounds for GPDs 
which have been obtained in [19]. 
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